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ABSTRACT 

To explain their observed radii, we present theoretical radius-age trajectories for the extrasolar giant planets 
(EGPs) TrES-4, XO-3b, and HAT-P-lb. We factor in variations in atmospheric opacity, the presence of an inner 
heavy-element core, and possible heating due to orbital tidal dissipation. A small, yet non-zero, degree of core 
heating is needed to explain the observed radius of TrES-4, unless its atmospheric opacity is significantly larger 
than a value equivalent to that at 10 x solar metallicity with equilibrium molecular abundances. This heating 
rate is reasonable, and corresponds for an energy dissipation parameter {Q p ) of ~10 3 8 to an eccentricity of 
^0.01, assuming 3xsolar atmospheric opacity and a heavy-element core of M c = 30 M$. For XO-3b, which 
has an observed orbital eccentricity of 0.26, we show that tidal heating needs to be taken into account to 
explain its observed radius. Furthermore, we reexamine the core mass needed for HAT-P-lb in light of new 
measurements and find that it now generally follows the correlation between stellar metallicity and core mass 
suggested recently. Given various core heating rates, theoretical grids and fitting formulae for a giant planet's 
equilibrium radius and equilibration timescale are provided for planet masses M p = 0.5, 1.0, and 1.5 Mj with 
a = 0.02-0.06 AU, orbiting a G2V star. When the equilibration timescale is much shorter than that of tidal 
heating variation, the "effective age" of the planet i s shortened, re s ulting i n evolutionary trajectories more like 
those of younger EGPs. Motivated by the work of iJackson et al.l d2008alfbh . we suggest that this effect could 
indeed be important in better explaining some observed transit radii. 

Subject headings: planetary systems — planets and satellites: general — stars: individual (HAT-P-1, GSC 
02620-00648, XO-3) 



1. INTRODUCTION 

As of the writing of this paper, an astounding 47 transit- 
ing extrasolar giant planets (EGPs) have been discovered 3 . 
For transiting planets, the inclination/planet-mass degeneracy 
is resolved and the photometric dip in the stellar flux dur- 
ing the transit yields the planet's radius (R n ). Theory then 
attempts to explain the measured radii dGuillot et al.l 19961: 



Burrows et all 12000. 2003, 2004, 2007; Bodenheim er et all 
200 11 120031: iBaraffe et all 120031 |200H 120081: IChabrier et al l 
2004lFortnev et al.ll2007HLaughlin et al.ll2005b . Importantly, 
the comparison between theory and measurement must be 
done for a given stellar type, orbital distance, planet mass, 
and age. The latter is poorly measured, but crucially impor- 
ta nt (see, e.g. g5 ]r e HP 209458b). 

iBurrows et all (12007b modeled the theoretical evolution 
of the radii of the 14 transiting EGPs known at the time. 
These authors suggest that there are two radius anoma- 
lies in the transiting EGP family, of which the smaller- 
radius anomaly is a result of the p r esence of dense core s 
(iMizunol [19801 iPollack et al.l 119961 iHubickvi et all I2004T) . 
whereas the larger-radius anomaly might be explained by 
the enhanced atmospheric opacities which slow down the 
heat loss of the core. They also discussed the effects on 
planet structure of the possible extra heat source in the in- 
terior, yet found no obvious correlation between the req- 
uisite power and the intercepted stellar power. Note that 
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3 See J. S chneider's Extrasolar Planet Encyclopaedia at 
http://exoplanet.eu the Geneva Search Programme at http://exoplanets.eu 
and the Carnegie/California compilation at http://exoplanets.org 



none of the transiting EGPs modeled by Bur rows et al.l (120071) 
is known to have a highly eccentric orbit. However, re- 
cently several transiting EGPs with significantly non-zero 
eccentricities (e > 0.15) have bee n discovered, including 
the EGPs XO- 3b Oohns-KruU et al.ll2008l). HAT- P-2b (a.k.a. 
HD 147 506b: iBakos et al.1 l2007at lLoeillet et al.1 120081). HP 
17156 b (Barbie ri et al.ll2007l: iGillon et al l 120081: iNarita et all 
120081: llrwin et alj l2008h7~and the "hot Neptune" GJ 436 b 
( Gillo n et al E007l:lDeming et al.ll2007l: iDemorv et al. l l2007l) . 

At least for those systems with hi ghly eccentric orbits, heat- 
ing d ue to orbital tidal dissipation (Bo denheimer et al]|2001l 
2003) should be incorporated into theoretical models for the 
radius-age trajectories. A preliminary exploration of this in a 
restricted context motivates the present paper. When the theo- 
retical radius evolution calculations are tailored to a system's 
specific planet mass, age, primary stellar properties, and or- 
bital distance (as they must), the current radii of most of the 
known transitin g EGPs can be e xpla ined by the theoretical ra- 
dius models of Bu rrows etail (120071) . In many instances, a 
higher-density core mass provides an even better fit, and ex- 
tra internal heat sources are not required. However, at least 
three transiting EGPs seem to be exceptional in some way. 
The planets TrES-4 and XO-3b (for its large-radius solution 
based on stellar parameters from s pectral synthesis model- 
ing) are cited by their discoverers (Mandushev et al.l 120071: 
iJohns-Krull et al.1 120081) as anomalously large and inconsis- 
tent w ith extant theoretical models. In addition. lBurrows et al.1 
d2007l) found that HAT-P-lb deviated from the core-mass 
stellar-metallicity relationship followed by many of the tran- 
siting EGPs that they modele d (for the core-m ass stellar- 
metallicity correlation, also see Guillot et al. 2001). 

Therefore, with this paper, we focus on this small subset 
of three interesting objects to determine how their radii can 
indeed be explained with minimal assumptions, that never- 
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theless can include tidal heating. We find that tidal heat- 
ing in the core, given measured (XO-3b) and possibly non- 
zero (TrES-4 and HAT-P-lb) eccentricities, which never- 
theless are still consistent with the upper limits, can nat- 
urally explain the measured radii. For HAT-P-lb with its 
new age estimate, we find that a core of reasonable size 
can now be accommodated. Importantly, when tidal heat- 
ing needs to b e invoked, a canonical tidal dissipation pa- 
rameter, "0„"(Gold reich & Soterl IT966) 4 . with a value near 
10 4-6 , along with the measured eccentricity (or reasonable 
values consistent with its current bounds), suffices to explain 
the measurements. Hence, only simple extensions of the de- 
fault evolution models that incorporate a known process with 
canonical parameters are required. We postulate that all mea- 
sured transiting EGP radii can be explained by available the- 
ory when proper account is taken of the measured planet- 
star system parameters, reasonable core masses that follow 
the relationship between core mass an d stellar metallicity 
dBurrows et al.li2007h iGuillot et al.ll2006l) . and tidal heating at 
the expected rate for non-zero, but measured, eccentricities. 
Exceptions to this might arise if it is determined that tidal dis- 
sipation occurs predominantly in the atmosphere, not the con- 
vective core, and/or if the eccentricity and semi-major axis 
history must be factored into the tidal heating history of the 
pl anet. The latte r effect is intriguing and has been suggested 
bv I Jackson etail (l2008allbb . 

We describe our computational techniques and model as- 
sumptions in ^2] In ^3] we review the measurements of these 
three EGP systems, identify the discrepancies between the ob- 
served planetary radii and those predicted by previous the- 
oretical models and present new theoretical radius-age tra- 
jectories for them using tailored atmospheric boundary con- 
ditions. These new trajectories and theoretical radii include 
the effects of tidal heating in the convective core for mea- 
sured or reasonable values of the orbital eccentricities and for 
a range of values for the tidal parameter, Q p . In $4] we pro- 
vide theoretical grids and fitting models for the equilibrium 
planetary radius and the equilibration timescale, given a cer- 
tain set of planet mass, orbital distance, and tidal heating rate 
for a G2V primary star. Finally, in <J3] we summarize our re- 
sults for each system, discuss the relevant constraints obtained 
on their structural properties, and list caveats concerning our 
model assumptions. 

2. COMPUTATIONAL TECHNIQUES AND MODEL ASSUMPTIONS 

A detaile d discussion of our computa tional techniques ca n 
be found in lBurrows et al.1 (120031) and lBurrows et al.1 (120071) . 
Here, we present only a brief summary, along with our model 
assumptions. We generate realistic atmospheres customized 
for the three EGPs, their time-averaged orbital separations, 
and primary stars. The adopted atmospheric boundary con- 
ditions incorporate irradiation using the observed stellar lu- 
minosity and spectrum, and the measured planet orbital dis- 
tance. For planets with eccentric orbits, the time-averaged 
insolation flux is employed in constructing the atmos pheric 
bound ary conditions. The theoretical stellar spectra of lKuruczl 
(1994) are adopted. For the given stellar spectrum and flux 
(inferred from the luminosity and the planet orbital distance), 
an S-T e ff-g grid is calculated for the core entropy, S, effec- 
tive temperature, T e ff, and gravity, g, using the discontinuous 
finite element (DFE) variant of the spectral code TLUSTY 

* Note that the Q„ p arameter used throughout this paper corresponds to the 
Q' p in lMardlind J2007l) . 



dHubenv & Lanzj[l995l) . It is assumed that both the stellar 
spectrum and flux are constant during the evol ution. 

We employ the Henyey evolutionary code of B urrows et al.l 
( 1 1997b for the radius-age evolutionary calculations, using the 
function T e ff(S,g) for the interior flux, inverted from the ta- 
ble of S, T e ff, and g, referred to above. The helium fraction 
Tee is assumed to be 0.25. We calculate models with differ- 
ent atmospheric opacities, the effect of which can be conve- 
niently mimicked by using lx solar, 3 x solar, and 10 x solar 
abundance atmospheres. Note that the increase in atmo- 
spheric opacity does not need to, and should not be due 
solely to, increased metallicity. The effects of increased atmo- 
spheric opacity and increased envelope heavy-element abun- 
dances are decoupled, so that the implied increases in the 
heavy-element burden of the envelope, if any, will not can- 
cel the expansion ef fect of enhanced atmospheric opacity 
dBurrows et al.| [2007). To model the presence of a heavy- 
element core, a compressible ball of olivine is placed in 
the center of the model planet, and pressure continuity be- 
tween the heavy-element core and the gaseous envelope is en- 
sured throughout the evolution. We adopt the Saumo net al.l 
(1995) equation of state (EOS) for the H2/ He envelope and 
the ANEOS by [Thompson & Lausor](ll972l) for olivine. 

As noted, to model the atmospheric opacity, we use su- 
persolar metallicities (e.g., 3 x solar, 10 x solar) to mimic 
the expansion effec ts of enhanced atmospheric opacity 
dBurrows et al.1 [2007). Possible causes for such enhanced 
opacities might be supersolar metallicities in the atmosphere, 
nonequilibrium chemistry, errors in the default opacities, and 
thick hazes or absorbing clouds. Note that the expansion ef- 
fects of enhanced atmospheric opacity and the shrinkage ef- 
fects of increased envelope metallicities conceptually decou- 
pled in our models and that an increase of the envelope heavy- 
element burden, will not necessarily cancel the expansion ef- 
fect due to enhanced atmospheric opacity. The effect of a 
central heavy-element core on the planet radius is to shrink it 
monotonically with core mass. 

Given a non-zero orbital eccentricity, the tidal dissipation 
rate is c alculated using the formulation summarized in th e Ap- 
pendix dBodenheimer et al.ll200lL 12003b iGu et al.ll2004T) . We 
assume that the planet is in synchronous rotation and that all 
the tidal heating is in the convective core. Note that the ef- 
fects of other core energy dissipation mechanisms on planet 
structural evolution are also implicitly addressed. As indi- 
cated in the Appendix, the tidal heating rate is proportional 
to /(e) / Q p , where /(e) = e 2 when e <C 1 . Also, the values 
of Q p for EGPs with masses M p ~ Mj, although very un- 
certain, are estimated to be ~ 1 4 — 1 7 (lAdams & Laughlin 
I2006t IGu et alJl2003t Hackson et alJl2008allbT) . Therefore, we 
calculate typical heating rates using the combination e 2 /Q p 
for TrES-4 and HAT-P-lb, of which the orbital eccentricities 
have been estimated to be <C 1, if nonzero at all. Models with- 
out external heat sources (e 2 /Q p = 0) are also presented. 

3. OBSERVED PROPERTIES AND THEORETICAL PLANETARY 
RADII 

Table Q] displays the relevant observed quantities and the 
corresponding references for the EGPs TrES-4, XO-3b, and 
HAT-P-lb and their parent stars. These properties include 
semi-major axis (a), orbital period and eccentricity (e), stel- 
lar mass (M*), radius (/?*), effective temperature T e ff, metal- 
licity ([Fe/H]*), age, and planetary mass (M p ) and radius 
(R p ). Also shown are the stellar flux at the planet's sub- 
stellar point, F p , in units of 10 9 erg cm" 2 s" 1 , and the ra- 
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tio between the possible tidal energy dissipation rate within 
the planet and the insolation rate £tide/2?insoiation in the unit 
of (gp/lO 5 )" 1 . Q- 1 = j^j /(-§)& is the specific dissipa- 
tion function of the planet, where Eq is the maximum energy 
store d in the tidal distortion a nd — ^ is the rate of dissipa- 
tion dGoldreich & Soteri 1 1 9661) The parameters for TrES-4 
and HAT-P-lb are drawn from lTorres et al.l ([2008). These au- 
thors provide a uniform analysis of transit light curves and 
stellar parameters based on stellar evolution models, and a 
critical examination of the corresponding error s. Since XO- 
3b has not been studied by iTorres et"afl (l2008|). we generate 
mo dels using paramete rs reported by Johns- Krull et al.l d2008l) 
and lWinn et all d2008l) . 

Theoretical evolutionary trajectories are presented for 
TrES-4, XO-3b, and HAT-P-lb under various assumptions 
about the atmospheric opacity, the presence of a heavy- 
eleme nt core and possible tidal heating. The transit radius 
effect dBurrows et al.ll2007h is also included in the models. 

3.1. TrES-4 

TrES- 4 is the current record - holder for the low est EGP 
density dMandushev et all 120071: ITorres et alj 120081) . In the 
discovery work, iMandushev et al .1 (120071) carried out spec- 
troscopic observations with the CfA Digital Speedometer 
(Lathaml 1 19921) , radial velocity (RV) measurements with 
Keck, and transit photometry in the z band with Kepler- 
Cam at the F. L. Whipple Observatory (FLWO) and in the B 
band using NASACam on the 0.8-m telescope at the Low- 
ell Observatory. Assuming [Fe/H]„ = O.O^qI dex, Mandu- 



shev et al. derived T f 
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1 73 8+0.092 p 
-1+0.10 



e/f 



6100!}^ K, 



1.223$ M Q , 



and an age of 4.7^ o Gy r f° r me star ' 



and M p = 0.84^;!" Mj and R p = 1 .674+, m 



+!ao94 Rj for the planet. 
Th e orbital eccentric i ty w as assumed to be exactly zero in the 
fit. Mandushev et al. (2007 ) suggested that its observed radius 
is too large to be exp lained by the theoretical EGP models of 
iBurrows et al.1 ([2007), given its estimated mass, age and inso- 
lation, even when the effects of higher atmospheric opacities 
and the transit radius correction are considered. 

The parameters of TrES-4 listed in Table [TJ are from 
ITorres et"ai . (2008). These authors derived these parame- 
ters using the RV measu rements and transit photometry from 
Mandushe v et al.1 (120071) . and the stellar atmospheric proper- 
ties from Sozzetti et al. (2008; in preparation). Note that their 
estimated planetary radius R p = 1 .75 l^o'oet * s larger 
than t he value of R p = 1.674^ °^ j; , f IMandushev et all 
d2007l) . 

3.1.1. Results for TrES-4 

We calculate the radius-age trajector ies for TrES-4 using 
the parameters from Tor res et al] ([2008), taking into account 
the possible effects of enhanced atmospheric opacities, and 
the presence of tidal dissipation given a small, yet non-zero, 
orbital eccentricity. Models with the presence of a heavy- 
element core are also calculated. The value of Zstide/Einsoiation 
quoted for TrES-4 in TableQ]is calculate d assuming e = 0.01 , 
but using all the other parameters from ITorres et all (120081) 
in which a circular orbit is assumed. More and better tran- 
sit observations are needed to better constrain its true orbital 
eccentricity. 

The first panel of Fig. [TJ shows the theoretical radii R p (in 
units of Rj) as a function of age (in units of Gyr) for TrES-4, 
under various assumptions concerning the atmospheric opac- 
ity and the level of tidal dissipation, without a heavy-element 
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FIG. 1. — Theoretical planet radius R p (Rj) versus age (Gyr) for TrES- 
4. Also sho wn on both panels with error bars are the observed radius 
and age from ITorres et all 120081) . Var ious values are assumed concern- 
ing the atmospheric opacity, the presence of a heavy-element core, and the 
core heating due to tidal dissipation. Different colors correspond to dif- 
ferent tidal heating rates, which are proportional to e~/Q p when c C 1 
(see the Appendix for more details). Top: This panel demonstrates the ef- 
fect of enhanced atmospheric opacity under various heating powers. Mod- 
els assuming 3 X solar ( 10 X solar) atmospheric opacities are plotted as solid 
(dashed) curves, whereas the black, blue, and red curves correspond to 
(e/0.01) 2 (10 5 /G P ) = 0.0, 0.1, and 6.0, all without a heavy-element core (M c 
= 0). Bottom: The effect of the presence of a heavy-element core is illustrated, 
where the dashed (solid) curves denote models with M c = 30 (0) Mm . If TrES- 
4 follows the core-mass stellar-metallicity relation found by Burr ows et all 
120071) . then it should contain a heavy-element core with M c ~20 — 40 A%, 
given its stellar metallicity [Fe/H]* (+0.14^^ dex). Assuming 3 X solar at- 
mospheric opacity, the model with Q p = 10 3 8 , e ~0.01, and a heavy-element 
core of M c = 30 Mm (black-dashed curve) explains the observed radius well. 



core (M c = 0). Models assuming 3x (lOx) solar atmospheric 
opacities are shown as solid (dotted) curves. The black curves 
show the models without any heat sources, whereas the blue 
(red) curves depict those with a tidal heating rate assuming 
(e/0.01) 2 (10 5 /2 P ) = 0.1 (6.0). Als o shown with error bars 
are the observed radius and age from lTorres et al.l d2008l) (Ta- 
ble [D- 

It can be seen from Fig. [JJthat if the radius and age esti- 
mates of ITorres et al.l (l2008l) do not deviate much from their 
true values (within the uncertainties), our models without any 
heat source s produce radii which a re ~3 a too small. As- 
suming the Mandush ev et al.l (120071) parameters, the theoret- 
ical radii are still ~2 a too small. The discrepancy will be- 
come smaller for models with even higher atmospheric opac- 
ities. So it is concluded that either the atmospheric opacity 
of TrES-4 is unusually large (much higher than the equivalent 
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of a 10 x solar metallicity, equilibrium mixture), or there are 
extra heat sources in the core. The required heating power 
is very modest; the model with (e/0.01) 2 (10 5 /G/>) = 0.1 and 
10 x solar atmospheric opacity produces theoretic al radii con- 
sistent with the 1-er lower bound of R p from iTorres et al.1 
(2008). 

The models shown in the first panel of Fig. Q]do n °t include 
any heavy-element cores. If TrES-4 follows the co r e-mas s 
stellar-metallicity relation studied by iBurrows et al.l (J2007), 
then there should be a heavy element core with M c —20-40 
M®, given its stellar metallicity [Fe/H]* (+0.14^;^ dex). The 
presence of such a heavy-element core will shrink the plane- 
tary radius, and, therefore, would require a higher tidal heat- 
ing rate than do models without a core to explain the ob- 
served radius. This effect of including a heavy-element core 
in the models for TrES-4 is shown in the second panel of Fig. 
Q] Assuming 3 x solar atmospheric opacity, the model with 
Q p = 10 3 8 , e —0.01, and a heavy-element core of M c = 30 M ffi 
{black-dashed) explains the observed radius well. Within 1 
a uncertainties, the model with Q p = 10 4 , e —0.01, and a 
heavy-element core of M c = 30 M® (red-dashed curve) can 
also fit the observed radius. 

In summary, unless the atmospheric opacity of TrES-4 is 
unusually large, core heating is required to explain its ob- 
served radius. However, the required heating power is mod- 
est. A non-core model with (<?/0.01) 2 (10 5 /Gp) = 0.1 pro- 
duces radii consiste nt with the 1 -a lower boundary of R p from 
ITorres et al l <|2008), assuming 10 x solar atmospheric opac- 
ity. The required energy dissipation rates become larger for 
models with a heavy-element core, but are still reasonable. 
For instance, the model with Q p = 10 3 8 , e —0.01, and a 
heavy-element core of M c = 30 M® produces the observed 
radius well, assuming 3 x solar atmospheric opacity. To bet- 
ter constrain models of TrES-4, definitive measurements of, 
or stronger limits to, its orbital eccentricity are needed. 

3.2. XO-3b 

XO-3b has been observed to be supermassive and on 
an eccentric orbit (M„ 
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iJohns-Krull et al.l 120081) . The discoverers obtained transit 
light curves with relatively small 0.3-m telescopes, spectro- 
scopic observations using the 2.7-m Harlan J. Smith (HJS) 
telescope and the 1 1-m Hobby-Eberly Telescope (HET) , and 
RV measurements with the HJS tele scope. Based on theoreti- 
cal spectral models of the HJS data. IJohns-Krull et all (120081) 
derive T eff = 6429tfg K, [Fe/H]* = -0.177^;^ dex, and 
logg* = 3.95!o;j]g2 f or tne star - Combined with the RV mea- 
surements, they arrive atM* = 1.41+^|M Q andi?* =2.13^;|j 
R G for the star, and M p = 13.25^;g M, and R p = 1.95^ 
R j for the planet. These authors have commented that XO- 
3b is observed to b e so large that in all cases analyzed by 
iFortnev et al.l (12007), their models predict a much smaller ra- 
dius. However, due to the absence of a precise trigonometric 
parallax of XO-3, its distance is very uncertain. Assuming a 
smaller distance, and, hence, a reduced stell ar mass and radius 
than o btained using the isochrone method, John s-Krull et al.l 
(2008) found a best fit to their transit light curves with logg* = 
4.19, M* = 1.24M Q , and R* = 1.48 R Q , with the correspond- 
ing estimates for M p of 12.03^ Mj and R p of 1 25$g Rj. 
These li ght-curve-based re sults have recently been strength- 
ened by IWinn et ail (120081) using larger aperture telescopes. 
These authors observed 13 transits photometrically using the 
1.2-m telescope at the FLWO, along with 0.4-0. 6-m tele- 
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FIG. 2. — Theoretical planet radius Rp (Rj) versus age (Gyr) for XO- 
3b. In the first panel, a planet mass of M„ = 13 .24 Mj based on the spec- 
tral synthesis model of Johns-Krull et al. (2008) and 3 X solar atmospheric 
opacity are assumed, whereas in the second panel the correspondin g values 
are M p = 11.79 Mj according to the light-curve fitting results of Winn et all 
12008) and 1 X solar atmospheric opacity. Also sho wn with error bars on both 
pan els are the observed radii and age estimates from Johns-Krull et al] <2008h 
and Winn et al. 12008), based on the two different analyses. In both pan- 
els, models assuming various tidal heating rates corresponding to Q p = oo 
(no heating), 1 7 , 10 6 , 10 5J and 1 5 5 are color coded. In all the mod- 
els, e = 0.260 ( Johns-Krull et al. 2008) is assumed. Note that the observed 
radii, based on the two different methods, differ quite a bit from one another 
(~50%). The very large radius based on spectral synthesis modeling can be 
fit by the model with Q p ~10 5,7 , whereas the much smaller one inferred from 
the light-curve fit can be expl ained by models with Q p values down to 10 6 , 
within l-cr errors. See 33.2. ll for more discussion. 

scopes. Based on these more precise transit light curves, they 
concluded that logg* = 4.244, M* = 1.213M , and R« = 1.377 
Rq, with the corresponding estimates for M p of 1 1 .79^ 59 Mj 
and for R p of 1.2 17^073 Rj. Since a trigonometric parallax 
measurement of XO-3 is still lacking that could distinguish 
these two different methods, we make models for both the 
spectroscopically-determined and light-curve-based parame- 
ter sets. 

Both the spectroscopical results of 1 Johns-Krull et aTl (T2 008) 
and the light-curve-based results of IWinn et alj d2008l) are 
listed in Table Q] We note that the M p and R p values derived 
using the two different methods differ significantly from one 
another (—10% in M p and —50% in R p ). Therefore, separate 
models and discussions for these different sets of planetary 
properties are in order and we calculate theoretical radii for 
XO-3b for both estimates of the planetary radius and mass. 

3.2.1. Results for XO-3b 
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We include the possible heating due to orbital tidal dissipa- 
tion, assuming reasonable values of Q p . Our results for XO- 
3b are shown in Fig. [2] where in the first panel a planet mass 
of M B = 13.24 Mj, derived from the spectral synthesis method 
of lJohns-Krull e t al. (2008), and an atmospheric opacity asso- 
ciated with 3 x solar metallicity are assumed, whereas in the 
second panel the corresponding value s are M„ = 1 1 .79 Mj 
from the light-curve fitting results of I Winn et al.l (120081) and 
1 x solar atmospheric opacity. In both cases, we assume that a 
massive heavy-element core is absent. For such massive plan- 
ets, the effect of reasonable core masses on planetary radii is 
small. 

As show n in the first panel of Fi g. |2] the very large radius 
derived bv lJohns-Krull et all (l2008). based on spectral synthe- 
sis modeling, can be explained by energy dissipation due to 
tidal heating when the Q p parameter of XO-3b is ~10 5 1 . For 
the case of the light-curve-based planetary mass (second panel 
of Fig. |2), the inferred radius is broadly consistent to within 
1-cr errors with the absence of an internal heat source, or with 
tidal heating with Q p > 10 6 . These specific constraints on Q p 
assume 3 x solar in the former case and 1 x solar atmospheric 
opacity in the latter. For our first model of XO-3b assuming 
M p =13 .24 Mj and 3 x solar atmospheric opacity, it would be 
surprising to find a planet with R p /Rj < 1.2 or > 2.2, since 
this will require Q p to be either too large or too small. For 
the second model assuming M p = 11.79 Mj and lx solar at- 
mospheric opacity, such a region would be R p /Rj < 1.1 or 
> 1.6, for the same reason. Given significant eccentricity, it 
is important to appropriately account for tidal heating in order 
to model the planet's structural evolution. We want to empha- 
size our adopted model assumptions that 1) the planet is in 
synchronous rotation and that 2) all the tidal heating is in the 
convective core. Even though the e 2 / Q p degeneracy is broken 
due to the known value of <?, detailed radius evolution models 
could be used to constrain Q p for EGPs, but only if it is deter- 
mined that tidal dissipation occurs predominantly either in the 
convective core or in the atmosphere, and if the uncertainties 
in the core mass and atmospheric opacity are both resolved. 

3.3. HAT-P-lb 

Using photometry conducted by the Hungaria n-made Au- 
tomated Telescope Network (HATNet) project, Bak os et all 
(2007b) discovered HAT-P-lb transiting one member of the 
stellar binary ADS 16402. These authors suggested that HAT- 
P-lb was too large to be explained by theoretical EGP models. 
Spectral synthesis modeling of the parent star ADS 16402B, 
based on its Keck spectra, yielded T e ff = 5975^ K and 
[Fe/H]„ = +0.13!g;jg dex. iBakos et all (l2007bl) also fit both 
stellar members in the binary to evolutionary tracks and based 
on the Subaru and the Keck spectra derived M» = 1 . 12±q '%Mq 
and R* = IA5^R Q for ADS 16402B and a best-fit age 
of 3.6 Gyr for th e binary. Using the z-band transit curves 
from KeplerCam dHolman et al.l 120061) . combined with RV 
measurements from Subaru and Keck, these authors derive 
M p = 0.53![J f 4 Mj and R p = 1 j^Rj, where the errors in the 
planetary radius include both statistical and systematic errors 
in both the stellar radius and mass. Note that in the above fits, 
a circular orbit ( e = 0) was assumed. However, the x 2 fitting 
of the RV data in Bakos et al. (2007b]) favors a small, yet non- 
zero, eccentricity: e = 0.09^q q|. These authors did estimate 
the heating rate assuming e = 0.09 and suggested that if this 
non-zero orbital eccentricity is confirmed, the observed large 
R p could be explained by tidal heating. However, note that 



a non-zero eccentricity is only suggestive, mainly due to the 
small number of RV observations (13 velocities, for which 
the typical S/N is about 150 per pixel). Since RV-based ec- 
centricity estimates are positively biased due to noise (e.g., 
IShen & Turnerl2008T), it is very likely th at the true e is smaller 
than 0.09. In fact, iJohnson et al.l (|2008) find a upper limit on 
e of 0.067, with 99% confidence, by combining their new and 
previous RV measurements. Therefore, we assume a smaller 
value, e = 0.01, in our baseline model and see where such an 
assumption leads. Such a small eccent ricity could result from 
Koza i cycles with tidal friction (e.g., iFabrvckv & Tremaind 
120071) . though t here is evidence th at the spin-orbit misalign- 
ment is small (IJohnson et al. 2008). 

B ased on more hig h-precision transit observations, how- 
ever, |WInnitil] (l20O7|) report that HAT-P-lb is less "bloated" 
than originally thought. Their observations include three tran- 
sits observed in z band with the 1 .2-m telescope at the FLWO, 
three observed through the "Gunn Z" filter (Pinfi ekTet al.l 
[19971) using the Nickel 1-m telescope at Lick Observatory, 
and three observed through the Johns on / filter u s ing th e 1- 
m telescope at the Wise Observatory. IWinn et all (12007b de- 

1 /3 

rived R*/Mj by fitting the transit light curves, and con- 
cluded that R* = 1.115^;^ Q and R p = 1 .203^ f 5 \Rj . Note 
that in their fits the orbital eccentricity was assumed to be 
zero. These authors suggest that the updated radius can b e 
explained by the structural models of iBurrows et al.l d2007l) . 
unless t he planet has a very m assive core of heavy elements. 
Indeed, IBurrows et al.1 (120071) calculated radius-age trajec- 
tories for HAT-P-lb. They included different core masses 
and atmospheric opacities in their models, and found that 
in order to fit the observed radius, HAT-P-lb deviates from 
the stellar-metallicity versus core-mass sequence otherwise 
roughly followed by the transiting EGPs included in their pa- 
per (see Fig. 9 of IBurrows et al.l 120071) . However, the stel- 
lar and planetary parame ters of the HAT-P-1 system adopted 
by IBurrows et alj ( 2007) are from the discovery work of 
iBakos et al.1 d2007bl). The p a ramete rs of the HAT-P-1 sys- 
tem derived by iTorres et al.1 (120081) are listed i n Table [D 
These authors compile the z-band light curves of I Winn et al.1 
d2007l) and the RV data and the atmospheric parameters of 
Bakos et al.l d200 7b). but with increased u ncertai nties for T e ff, 
[Fe/H]*, and logg*. As we show in 33.3.11 using these 
new parameters and the new planet radius, we now find in- 
ferred core masses that are roughly consistent with the stellar- 
metallicity vers us core-mass relation ship followed by the 
EGPs studied by IBurrows et all (120071) . 

3.3.1. Results for HAT-P-lb 

Give n the new measured radii and stellar age of lTorres et al.l 
(2008), we have reexamined the best-fit core masses for HAT- 
P-lb. The effects of the possible heating due to tidal dissi- 
pation on the planet's structural evolution are considered, as- 
suming a small yet non-zero eccentricity e = 0.01 . Theoretical 
evolutionary trajectories of planet radius with age for HAT-P- 
lb are shown in Fig. [3] where the first (second) panel displays 
the results assuming a 3 x solar ( 10 x solar) atmospheric opac- 
ity. Models with different tidal heating rates proportional to 
(e/0.01) 2 (10 5 /2 P ) are color-coded as labeled. For both of 
the panels, different line styles represent "no heavy-element 
cores," or the presence of a heavy-element core with a range 
of masses in units of Earth masses, M ffi . For clarity, only se- 
lected models are presented in the figure. Table|2]lists the best 
estimates for the core mass under various assumptions. 
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FIG. 3. — Theoretical planet radius R p (Rj) versus age (Gyr) for HAT-P-lb. 
Also sho wn with error bars are the observational stellar age and planet ra- 
dius from Torres et al. (2008), as listed in Table[TJ We reexamine the best-fit 
core masses for HAT-P-lb, under various assumptions concerning the atmo- 
spheric opacity and the possible core heating rate due to tidal dissipation. 
The top (bottom) panel shows models assuming 3 X solar ( 10 X solar) atmo- 
spheric opacities. In both panels, the black, blue, and red curves denote 
models assuming 

(e/0.01) 2 (10 5 /G P ) = 0, 0.8, and 8, respectively. Differ- 
ent line styles correspond to models with various heavy-element core masses 
in units of the Earth mass, as labeled on the plot. If there is tidal heat- 
ing assuming reasonable values of (e/0.01) 2 (10 5 /2p) ~0.8, then the core 
mass required to fit the observed radius is ~30 M9 (~40 Mg), assuming 
3 X solar ( 10 X solar) atmospheric opacity, in which ca se HAT-P-lb does fol - 
low t he cor e-mass stellar-metaliicity relation found by Burrows et al] 120071) . 
See 33.3. ll for a discussion. 

As demonstrated in Fig. [3] and Table [2] there are multiple 
solutions to explain the observed radius of HAT-P-lb. In all 
the cases considered, a non-zero heavy-element core mass is 
needed, which, without any external heat sources, is —15 M ffi 
(—20 M®) assuming 3xsolar (lOxsolar) atmospheric opac- 
ity,. This best-fit core mass becomes larger when there is 
external heating. If HAT-P-lb does follow th e approximate 
core-m ass/stellar-metallicity relation found by iBurrows et al.l 
(2007), then given its [Fe/H]„ (+0.13!oog dex), the core mass 
would be —30 M e (—40 M ffi ) assuming 3xsolar (lOxsolar) 
atmospheric opacity. These core masses correspond to the 
cases with (e/0.01) 2 (10 5 /2/>) = 0.8. 

In summary, if there is tidal heating with reasonable val- 
ues of Q p , the core-mass estimates suggest that HAT-P-lb fol- 
lows the correlation between s tellar metallicity and core mass 
found by IBurrows et al.l (|2007), or if there is no extra heating, 
deviates mildly from the correlation sequence. However, a 
lar ger core mass, more i n keeping with the correlation found 
by IBurrows et al.l (120071) . is more consistent with reasonable 



values of (e/0.01) 2 (10 5 /2/:>X as long as e is non-zero and Q p 
is not anomalously small. 

4. EQUILIBRIUM PLANETARY RADII AND EQUILIBRATION 
TIMESCALES FOR VARIOUS HEATING RATES: A PARAMETER 
STUDY 

In this section, we investigate the effects of generic core 
heating on the planet's equilibrium radius, and the time to 
reach this equilibrium. We refer to the latter as the equilibra- 
tion timescale. The heat source discussed in this section could 
be due to orbital tidal heating, but does not have to be. Figure 
H shows equilibrium planetary radii (R eq ) and equilibration 
timescales (j eq ) with various core-heating powers, quantified 
by Ei ieating /E imo i at i orl . Models are calculated for planets with 
masses M p = 0.5, 1.0, and 1.5 My, with semi-major axes a = 
0.02, 0.03, 0.04, 0.05, and 0.06 AU, orbiting a G2V star. Equi- 
librium is defined as the state after which the planet radius is 
constant to within a part in 10 5 . We define the equilibration 
timescale as the time it takes the planet to evolve from 1.25 
R eq to 1.05 R eq . Note that for those rare models for which the 
planetary radii still change by more than 10~ 5 R eq at the end of 
calculation (10 Gyr), the equilibrium state is assumed to have 
been reached at the final age of the evolutionary trajectory. 

In Fig. |H filled circles represent the values of R eq 
and r eq calculated from our theoretical trajectories, whereas 
the curves are least-square fits to them. We adopt a 
fourth-order polynomial in fitting R eq /Rj as a function of 
\og(E healin g/E imnlarion ), and a linear model for log(r e? /Gyr) 
versus \og{E hea , ing /E imolation ), given by: 

^ = C + Cix + C 2 x 2 + C3X 3 + C 4 x 4 , 

( r \ (1) 

( = b + k r 



log 



\Gyr/ 



1 .25 Rj and T eq — 2 Gyr for 
heating! ^insolation ^; ^ LlJ « e , ~ 1 -8 Rj and T eq — 3 Myr for 



where x = \og(E„ eating /Ei„ so i a , io „). The model fits of the pa- 
rameters are given in Table [3] 

For an extreme close-in EGP with M p = 0.5 Mj at a = 0.02 
AU orbiting a G2V star, the equilibrium planetary radii range 
from ~1.3 Rj for little heating {E heating /E inso , ation < 10" 6 ) to 
-2.5 Rj for strong heating (E heating /E insolatio „ ~ 10" 3 ). The 
corresponding timescales for the planet to reach these radii 
are —3 Gyr for the former and —3 Myr for the latter. At 
a = 0.06 AU, the equilibrium radii are smaller and the relevant 
timescales are longer - from R eq 

E Stealing I E Insolation ^ 10 tO R. 

E heating/ 'E 'insolation ~ 10~ 2 . For more massive planets, the equi- 
librium radii are smaller and the timescales are longer. For an 
EGP with M p = 1.5 Mj at a = 0.02 AU, the values are R eq ~ 1 .8 

Rj and r eq - 5 Myr for E hea ting/Ein S oiation ~ 10" 2 . Our theoret- 
ical model grids along with the fitting curves provided in eq. 
(Q~|) and the parameters listed in Table [3] can be used to calcu- 
late the equilibrium planetary radius and the typical timescale 
to reach it, given different combinations of planet mass, or- 
bital distance, and the ratio of core-heating power to insola- 
tion power. 

5. SUMMARY AND DISCUSSION 

We have calculated theoretical radius-age trajectories for 
three EGPs: TrES-4, XO-3b, and HAT-P-lb, under various 
assumptions concerning atmospheric opacity, the presence of 
an inner heavy-element core, and possible heating due to or- 
bital tidal dissipation. The main model results are the follow- 
ing: 
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FIG. 4. — Equilibrium planetary radius R eq (Rj) and equilibration timescale (Gyr) assuming various ratios between the core heating power and the insolation 
power. The equilibration timescale is defined as the time it takes the planet to evolve from 1.25 R eq to 1.05 R eq . Models are calculated for planets with masses 
M p = 0.5, 1.0, and 1.5 Mj, and semi-major axes a = 0.02, 0.03, 0.04, 0.05, and 0.06 AU, orbiting a G2V star. Filled circles represent results calculated from 
radius-age trajectories, whereas the curves are fits to them given by eq. (TJ and the corresponding parameters in Table[5] See ^4]for more information. 

no core heating sources, or with tidal heating assum- 
ing Q p > 10 6 , within 1-cr errors. These constraints on 
Q p assume 3 x solar atmospheric opacity for the former 
case and 1 x solar for the latter, but are only weakly de- 
pendent on this. 

3. We have reexamined the core mass required for HAT-P- 
lb using the update d data (importantly, its radius) from 
iTorres et"aT] (|2008), and now find it generally follows 
the correlation between core mass a nd stellar metal- 
licity found by iBurrows et alj d2007t) . In all the cases 
considered, a non-zero heavy-element core mass is 
nee ded to explain the observed radius (R p = 1 .24-2^ q|| 
Ri: ITorres etail 120081) . The core mass is ^15 M ffi 
(~20 M ffi ) assuming 3xsolar (lOxsolar) atmospheric 
opacity when there is no external heating. If there 
is tidal heating corresponding to reasonable values of 
(e/0.01) 2 (10 5 /2 p ) ~0.8, then the core mass required to 
fit the observed radius is ~30 M ffi (~40 M ffi ) assuming 
3 x solar (lOxsolar) atmospheric opacity, in which case 
HAT-P-lb follo ws the core-mass stel lar- metallicity re- 
lation found by IBurrows et al l (120071) and lGuillot et all 
d2006h . 

In addition, we have carried out a parameter study of the ef- 
fects of core heating and provided theoretical grids and fitting 
formulae for the equilibrium planet radius and equilibration 
timescale, given various core heating powers for planets with 
masses M p = 0.5, 1.0, and 1.5 Mj with a = 0.02-0.06 AU, 
orbiting a G2V star. The fitting formula for the equilibrium 
planet radius can be used for a theoretical zeroth-order esti- 



1. Unless the atmospheric opacity of TrES-4 is unusu- 
ally large (much higher than 10 x solar equivalent), core 
heating is required to explain its observed radius {R p = 

L751 -o:o62 R f> ITorres et al.ll2008l) . However, the re- 
quired heating power is modest. A non-core model with 
(e/0.01) 2 (10 5 /£2p) = 0.1 produces radii consistent with 
the 1-cr lower boundary of R p from lTorres et alj ((2008), 
assuming lOxsolar atmospheric opacity. The required 
energy dissipation rates become larger for models with 
a heavy-element core. The model with Q p = 10 3 8 , 
e ^0.01, and a heavy-element core of M c = 30 M®, re- 
produces the observed radius well, assuming 3 x solar 
atmospheric opacity. If TrES-4 follows the core-mass 
stellar -metallicity correlation found by IBurrows et al.1 
(2007), then the models with a non-zero heavy-element 
core mass are favored, consi dering its stellar m etallic - 
ity [Fe/H]„ = +0.14^$ dex (ITorres et al.l2008h . Ongo- 
ing Spitzer photometry of its secondary eclipse will put 
more stringent constraints on e and can either confirm 
or rule out these possibilities. 



For XO-3b, we have shown that orbital tidal heat- 
ing is a key factor in explaining the planet radius. 
The very large radius (R „ = 1.95+[] Rj) derived by 



iJohns-Krull et alj d2008l) based on spectral synthesis 
modeling can be explained by energy dissipation due to 
tidal heating. In this case, the Q p parameter of XO-3b 
is near ^10 5 7 , a not unreasnoable value. On the other 
hand, the much sma ller radius (1.21 7^" n7 3 Rj) based on 
light-curve fit by IWinn et al] d2008l) is consistent with 



8 



LIU ET AL. 



mate, without carrying out detailed evolutionary calculations. 
The equilibration timescale r eq characterizes the time it takes 
the planet to adjust its structure in response to a given degree 
of core hea ting. 

Recently Jack son et alJ (l2008ah considered the effect of the 
co-evolution of the orbital eccentricity and the semi-major 
axis on the tidal dissipation history. In the past, the semi- 
major axis had been assume d to be constant when conducting 
tidal evolution st udies (e.g. iBoden heimer et al. 2001, 12003c 
lGuetalJl2004h . Uackson et alJ d2008bl) calculate the evolu- 
tionary histories of the tidal dissipation rate for several EGPs, 
and find that in most cases the tidal heating rate increases as 
a planet moves inward and then decreases as the orbit circu- 
larizes. The relevant timescale, r/ !efl „„ g , is the time it takes the 
tidal heating rate to decay by a factor of e. If r eq ^> Treating, 
then it is valid to take a constant effective tidal heating rate 
in the planet radius-age trajectory calculation. If r eq ^Th ea ting, 
then in order to account for the effect of a varying tidal heat- 
ing rate, different values should be used at each time step of 
the radius-age trajectory calculation. If r eq <C Treating, then the 
planet will have enough time to reach an equilibrium radius 
before the tidal heating rate decays significantly. In this case, 
the planet's structure evolves in a quasi-equilibrium manner. 
Theoretical planet radius-age trajectory models will easily be 
able to account for the effect of varying tidal heating rate by 
adopting different core heating rates at each time interval of 
Theating during the calculation. In effect, there is a reset of the 
"clock right after the time of maximum heating - the planet 
becomes most extended on a timescale ~r e9 after the tidal 
heating rate achieves this maximum. Because of the intense 
heating and the quick response, the planet loses the memory 
of its shrinkage history before maximum heating, which is ef- 
fe ctively a reset of its "a ge .' ' 

iBurrows et alJ (120071) calculated theoretical radii for HD 
209458b. They found that the measured radius deviated at 
the ~ 1.5-cr level for the age they assumed, even when em- 
ploying 10 x solar atmospheric opacity, no inner solid core, 
and no core heating. However, the updat ed age measure- 
ment for HD 209458b bv lTorres et all j2008l) of 3.1^ Gyr is 
much smaller than the one adopted by IBurrows et al.l (2007) 
(5.5^["5 Gyr), whereas the updated radius, 1. 359^0 0)9, is sim- 
ilar to that used by IBurrows et all d2007l) (1.32^). As a 
result, the iTorres et all Q2008) radius and age measurement 



for HD 20945 8 b can be fit by the 10 x solar-opacity model of 
IBurrows et alJ ((2007) within — l-o uncertainties, without the 
need of any core heating sources. Moreover, if the tidal heat- 
ing rate of HD 209458b decayed fro m -4xl0 26 erg s" 1 t o 
— 4x 10 24 erg s" 1 during the past 2 Gyr dJackson et alJ2 008b). 
then its effective age would be 2 Gyr younger, due to the 
"clock reset" effect. Based on our parameter study in $31 
the equilibration timescale of HD 209458b during maximum 
heating would have been —0.05 Gyr. This is small enough 
compared with the decay timescale of the tidal heating rate 
(-1 Gyr from — 4xl0 26 erg s" 1 to -1.5 xlO 26 erg s" 1 ) for 
our models t o reproduce HD 20945 8b 's current radius. In 
this case, the ITorres et al.l (120081) radius measurement for HD 
209458b, along with an "effective age" of —1.1 Gyr, can ex- 
plain HD 209458b's radius within — l-o uncerta i nties, even 
with the 1 x solar-opacity model of Bur rows et alJ (120071) and 
a core mass of 10-20M®, but without the need for any current 
core heating. The latter comports with the very small limit of 
—0.001 on its current orbital eccentricity. A p ossible caveat 
to the tidal evolution scenarios described in iJackson et alJ 
(2008a b) is that they can be dramatically cha nged due to eve n 
a small undetected perturbing body (e.g., Mar dlindl2007l) . 
Sensitive searches for such companion bodies are needed to 
further constrain this possibility. Nevertheless, future stud- 
ies should combine orbital semi-major axis, eccentricity, and 
planet radius evolution models in a more coupled fashion. It 
is not only more consistent to consider the co-evolution of a, 
e, and R p , but also important to factor in the feedback of the 
associ ated tidal heating po wer on R p and its radius-age trajec- 
tory (I Jackson et al.ll2008allbh . Such a project is in progress. 
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APPENDIX 



EXTERNAL HEATING DUE TO TIDAL DISSIPATION 



The total tidal energy dissipa tion rate within the planet in its rest frame assuming equilibrium tides with constant lag angle and 
synchronous rotation is (e.g. lGoldreich & Soterfll966t rBodenhei mer et al.|[200Tl l2003h iGu et al.ll2004 : 



^tide 



GM*UM. 

flfcirc 



1.1 x 10 24 ergs _I 



0.01 



Mp 
M, 



-1/ T . N-l 
/ ' circ 



0.05 AUy \Gyr 



(Al) 



where fi=M r M p /(M r +M p ) is the reduced mass, /(e) = ^[h^(e)-2h4(e)+h 5 (e)] is a function of orbital eccentricity with /13(e) = 
(1 + 3e 2 + 3e 4 /8)(l - e 2 )" 9 / 2 , h A (e) = (1 + 15e 2 /2 +45e 4 /8 + 5e 6 /16)(l - e 2 )" 6 , and /15(e) = (1 + 31e 2 /2 + 255e 4 /8 + 185e 6 /16 + 
25e 8 /64)(l -e 2 )" 15 / 2 dGu et al.ll2004l) . Note that /(e) — > e 2 as e — > 0. r c i rc denotes the circularization timescale, which is 
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A more informative quantity is the ratio of the tidal energy dissipation rate and the insolation rate, given by 

gtide _. GM tf j,f(e) 

^insolation ^^P^p aT ckc 

~ 69xrW— V[^l f Q " VY M * Y /2 ( R »Y( a V 15/ V Fp X 

VO.01/ [ e 2 J V 10 5 / \MqJ \RjJ V0.05AUy/ \ 10 9 ergcm- 2 s" 1 J 
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TABLE 1 

Observational Properties of the Transiting Planet Systems. 



a Period M, R, T eff [Fe/H]» Age M„ R p F p * E Me /E mmllltmn b 

System (AU) (days) e (Mq) (R©) (K) (dex) (Gyr) (Mj) (Rj) (10 9 erg cm" 2 s" 1 ) KQp/l&T'] References 

xn _-) nn47fi +00005 1 1 Q 1 54. n?fin +0017 ] 41 +0.08 2 13+0- 21 6429+ 50 -0 177+ 0027 2 69+ 014 n 25+ - 64 i 95+O.I6 4 20+0-84 SQylfr 1 c 

AU J u - u4/D -0.0005 J-iyiJt U.zou„ 017 i-^i-o.os Z ' 1J -0.21 O4zv -50 u - 1 "-0.027 zoy -0.I6 1J - z;, -0.64 1 - : "-0.16 4 ' zu -0.84 x,xlu 

n n454+ 000082 119152 n2fin +0017 1 2n +0066 1 <!77+0-083 fi429 +100 -0 177+0-080 2R2+ - 58 1179+0- 551 I 217+ - 073 1 91+0-27 5 1 y ifl-l d 
U - U4:>4 -0.00082 J-iyi3Z U.ZOU„ 017 l zl -3_0.066 1 " J/ '-0.083 O4Zy -100 u - 17 '-0.080 Z - BZ -0.82 11 ' 7y -0.59 '-0.073 1 - y; '-0.27 3.1 X 1U 

TrES-4 0.05092^ 00072 3 553945 „ .0 1.394^;^g 1.816f^;gg 6200^5 +0.14+O;0S 2.9+0; 4 0.920+0; , 73 ; 1.751^ 0g 2.31t8;IA 2.3 x 1(T* e f 

HAT-P-1... 0.0553;000i2 4.45543 _ 00 l.^+O; , 7 ^ 1.135+0; 048 5975+}2° +0.13! ,; ! 2.7! 2 ; ' O.532;0-030 1.242^ 0.66+O; 08 9.2 x lO" 5 <= f 



a The stellar flux at the planet's substellar point.k is the total tidal energy dissipation rate within the planet in its rest frame (Gu et al. 2004). insolation = wRjFp- iJoh ns-Krull et al. 1 2008); Inferred from spectroscopically derived stellar 
parameters.^ Winn et al. 1 2008); Determined from light-curve fits. e Assuming e = 0.01 / Torres et al. 1 2008); Note that they assume e = exactly in deriving the parameters, since the radial-velocity data are consistent with a circular orbit. 



THEORETICAL RADII OF EXTRASOLAR GIANT PLANETS 



TABLE 2 

Best-estimate Core Masses for HAT-P- 1b under 
Various Assumptions. 

Atmospheric Opacity 3 X solar 10 X solar 

, 2 



10? 

0.01 I \ Q. 



0.8 8 O.i 



Mc/M® 15 30 60 20 40 80 

NOTE. — The tidal dissipation rate is proportional to e 2 /Q p when 
e « 1. 



TABLE 3 

Fitting Parameters for the Equilibrium Planetary Radius and 
Equilibration Timescale for EGPs Undergoing Core Heating. 



M p a 



(Mj) 


(AU) 


Co a 


Ci a 


c 2 a 


c 3 a 


c 4 a 


b b 


k b 


0.5 


0.02 


15.4 


9.69 


2.63 


0.328 


0.0156 


-4.78 


-0.845 




0.03 


7.74 


4.43 


1.22 


0.157 


0.00785 


-4.26 


-0.794 




0.04 


6.42 


4.00 


1.25 


0.182 


0.0101 


-3.95 


-0.762 




0.05 


4.61 


2.66 


0.866 


0.135 


0.00809 


-3.77 


-0.746 




0.06 


3.60 


1.79 


0.573 


0.0905 


0.00560 


-3.71 


-0.758 




0.02 


6.42 


4.51 


1.58 


0.255 


0.0155 


-4.10 


-0.783 




0.03 


4.44 


3.02 


1.16 


0.208 


0.0139 


-3.70 


-0.742 




0.04 


3.42 


2.19 


0.894 


0.170 


0.0122 


-3.42 


-0.701 




0.05 


2.72 


1.54 


0.669 


0.137 


0.0106 


-3.27 


-0.686 




0.06 


2.33 


1.15 


0.512 


0.110 


0.00897 


-3.17 


-0.677 




0.02 


4.20 


2.54 


0.888 


0.145 


0.00899 


-3.59 


-0.702 




0.03 


3.15 


1.76 


0.663 


0.118 


0.00796 


-3.21 


-0.654 




0.04 


2.57 


1.34 


0.543 


0.105 


0.00762 


-3.01 


-0.630 




0.05 


2.13 


0.924 


0.391 


0.0808 


0.00634 


-2.90 


-0.622 




0.06 


1.91 


0.693 


0.292 


0.0610 


0.00486 


-2.82 


-0.612 



a R eq /Rj = C +C i x + C 1 x 1 +C i x i +Ca 4 , where x = log(E kealing /E ima , alim ). h log(-r e ,/Gyr) = 
b + kx. 



